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The  Earth's  gravitational  potential  is  now  usually  expressed  in  terms  of  a 
double  series  of  tesseral  harmonics  with  several  hundred  terms,  up  to  order  and 
degree  at  least  20.  The  harmonics  of  order  14  can  be  evaluated  by  analysing 
changes  in  satellite  orbits  which  experience  14th-order  resonance,  when  the 
track  over  the  Earth  repeats  after  14  revolutions. 

In  this  Report  we  describe  our  first  evaluation  of  individual  l4th-order 
coefficients  in  the  geopotential  from  analysis  of  the  variations  in  inclination 
and  eccentricity  of  satellite  orbits  passing  through  14th-order  resonance  under 
the  action  of  air  drag.  Using  results  from  eleven  satellites,  we  find  the 
following  values  for  normalized  coefficients  of  harmonics  of  order  14  and 


degree  l , C 


Sfc, 14  ’ for 


*.  - 14,  15  ....22: 


These  values  provide  a test  of  the  accuracy  of  the  14th-order  coefficients 
in  comprehensive  geoid  models.  Detailed  comparisons  with  three  recent  models  are 
made,  showing  good  agreement  on  some  coefficients  and  discrepancies  on  others. 
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1 INTRODUCTION 

One  of  Che  finest  achievements  in  the  Earth  sciences  in  the  past  few  years 
has  been  the  accurate  determination  of  the  shape  of  the  geoid,  obtained  by 
assuming  that  the  geopotential  can  be  expanded  in  a double  series  of  tesseral 
harmonics,  then  analysing  hundreds  of  thousands  of  satellite  observations  and 
solving  the  equations  to  obtain  values  for  hundreds  of  the  harmonic  coefficients. 
Examples  of  such  comprehensive  geoid  models  already  published  include  the 
Goddard  Earth  Model  8 (GEM  8)*,  the  Smithsonian  Standard  Earth  IV. 3 (SSE  IV. 3)^ 
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and  the  European  GRIM  2 model  . The  accuracy  of  the  solutions  has  been  improving 
rapidly  as  more  laser  ranging  data  are  included:  the  most  recent  Goddard  Earth 
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Model  , GEM  10,  includes  598  geopotential  coefficients  and  is  complete  to  degree 
and  order  21;  GEM  10  utilizes  840000  observations  including  200000  laser  ranges 
(many  on  Geos  3 and  Lageos),  and  gives  a geoid  which  may  be  accurate  to  about 
1 metre  over  most  of  the  world. 


Though  the  overall  geoid  is  well  determined  in  these  comprehensive  solu- 
tions, most  of  the  values  of  individual  geopotential  coefficients  of  degree 
higher  than  about  11  are  poorly  determined:  in  GEM  7,  Wagner  estimated  that 
the  errors  exceed  30Z  for  degrees  higher  than  12.  The  values  of  coefficients  of 
a particular  order  can,  however,  be  accurately  determined  by  analysing  the  orbits 
of  satellites  which  experience  resonance  of  that  order,  as  their  orbits  slowly 
contract  under  the  influence  of  air  drag.  For  example,  15th-order  resonance 
occurs  when  the  Earth  spins  once  relative  to  the  orbital  plane  while  the  satellite 
completes  15  revolutions.  The  satellite  then  crosses  the  equator  at  intervals  of 
24°  in  longitude,  and  its  track  over  the  Earth  repeats  after  15  revolutions.  The 
l5th-order  harmonics  in  the  geopotential  have  maxima  every  24°  in  longitude,  and 
therefore  have  the  same  effect  on  the  satellite  on  every  revolution.  Consequently 
the  perturbations  which  they  produce  tend  to  build  up,  and  if  near-resonance 
continues  for  several  months,  the  resulting  changes  in  the  orbital  parameters 
(particularly  the  orbital  inclination)  are  large  enough  to  be  analysed,  and  to 
give  an  accurate  value  for  a 'lumped  harmonic’  appropriate  for  that  orbit  - the 
lumped  harmonic  being  a linear  function  of  the  individual  harmonics  of  order  15 
and  degree  15,  17,  19  ....  (if  the  inclination  is  being  analysed).  By  obtaining 
such  values  from  orbits  at  differing  inclinations,  the  values  of  the  individual 
coefficients  can  be  determined.  The  15th-order  resonance  is  particularly  fruit- 
ful, because  it  occurs  when  the  average  height  of  a satellite  over  the  Earth  is 
about  500  km,  and  numerous  satellites  are  'dragged  through'  15th-order  resonance. 
Three  years  ago  we  analysed  the  changes  at  15th-order  resonance  for  13  satellites. 


obtaining^’ ^ values  for  coefficients  of  order  15  and  degree  15,  17,  19  ....33 
(from  variations  in  inclination)  and  coefficients  of  degree  16,  18,  20  and  22 
(from  variations  in  eccentricity).  These  sets  of  values  for  15th-order  harmonic 
coefficients  were  better  than  the  conflicting  sets  in  the  comprehensive  geoid 
models  existing  at  that  time,  and  have  subsequently  been  incorporated  in  GRIM  2. 

The  success  of  this  analysis,  which  also  yielded  tentative  values  for 
30th-order  harmonics,  encouraged  the  idea  of  studying  other  resonances.  We 
looked  at  the  feasibility  of  utilizing  'half  resonances',  such  as  29:2  (when  the 

g 

satellite  makes  29  revolutions  while  the  Earth  spins  twice),  and  in  1976  Walker 
succeeded  in  analysing  29:2  resonance  for  Ariel  1,  and  obtained  values  for  lumped 
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29th-order  harmonics;  while  Hiller  and  King-Hele  analysed  31:2  resonance  for 
Proton  4,  though  with  less  accurate  results,  because  of  the  high  drag. 

The  other  obvious  resonances  to  study  are  14th  and  16th  order.  The  latter 
is  difficult  because  16th-order  resonance  occurs  when  the  average  height  of  the 
satellite  is  near  200  km,  when  the  drag  is  usually  severe,  and  the  resonance  does 
not  last  long  enough  to  build  up  a measurable  perturbation.  For  14th-order 
resonance  the  average  height  is  800  km,  and  the  problem  is  lack  of  drag:  the 
satellites  linger  on  the  brink  of  resonance  for  many  years  before  experiencing 
the  effects  fully.  After  waiting  for  two  years,  we  now  have  four  satellites  which 
give  good  results  at  14th-order  resonance  and,  by  combining  these  with  three 
analysed  by  Wagner^  and  results  from  four  high-drag  satellites,  we  have  evaluated 
14th-order  coefficients  of  degree  14,  15,  16  ....22,  with  tentative  values  for 
higher  degrees. 

2 THEORETICAL  EQUATIONS  FOR  ORBITAL  CHANGES  NEAR  1 4TH-0RDER  RESONANCE 

If  we  accept  that  a double  infinite  series  of  tesseral  harmonics  is  the  best 
representation  of  the  geopotential, its  longitude-dependent  part  can  be  written 
in  normalized  form  at  an  exterior  point  (r,0,A)  as 

00  t, 

r Yj  Z (rf  (C0S  6)  l5*m  COS  mX  + 5*m  sin  “’Km  ’ (,) 

im2  m**l 

where  r is  the  distance  from  the  Earth's  centre,  6 is  co-latitude,  X is 

longitude  (positive  to  the  east),  p is  the  gravitational  constant  for  the 
3 2 

Earth  (398601  km  /s  ),  R is  the  Earth's  equatorial  radius  (6378.1  km), 

P™  (coa  6)  is  the  associated  Legendre  function  of  order  m and  degree  l , 
and  C and  S are  the  normalized  tesseral  harmonic  coefficients,  of  which 
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those  of  order  m ■ 14  particularly  concern  us  here.  The  normalizing  factor 

N is  given  by11 
im 


N 


Im 


2(21  + l)(t  - m)l 
(£  + m) ! 


(2) 


The  rate  of  change  of  inclination  i caused  by  a relevant  pair  of 

_ - j 2 

coefficients,  C.  and  S.  , near  B:a  resonance  may  be  written 
im 


im 
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where  F.  is  Allan's  normalized  inclination  function  , G.  is  a function 

v imn  * Dnn 


Imp 


£pq 


- _ j 2 

of  eccentricity  e for  which  explicit  forms  have  been  derived  by  Gooding  , A 


denotes  'real  part  of'  and  j ■ /-T  . The  resonance  angle  9 is  defined  by 
the  equation 

* - a((u  + M)  ♦ - v)  , (4) 


where  u is  the  argument  of  perigee,  M the  mean  anomaly,  fl  the  right 


ascension  of  the  node  and  v the  sidereal  angle.  The  indices  y,  q,  k and  p 


m equation  (3)  are  integers,  with  y taking  the  values  1,  2,  3 ....  and  q 
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the  values  0,  ±1,  ±2,  ....;  the  equations  linking  l,  m,  k and  p — 


are 


yB  ; k - ya  - q ; 2p  ■ l - k . 


1 , and  the  m-suffix  of  a relevant  (C.  ,S  ) pair 

Jem  £m 


Here  8 - 14  and  a 
is  given  uniquely  by  the  choice  of  y . The  values  of  l to  be  taken  must  be 
such  that  i % m and  (£  - k)  is  even.  The  successive  coefficients  which  arise 
and  q ) may  usefully  be  gathered  together  in  a lumped  form  and 


(for  given  y 
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written  as 
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where 
the  Q 


l l 

l increases  in  steps  of  2 from  its  minimum  permissible  value  l. 


and 


are  constant  coefficients  with  Q 
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For  the  14:1  resonance  the  most  important  terms  in  equation  (3)  are  likely 


to  be  those  with  y 


J 


because  y * 2 gives  m - yB  ■ 28  , and  the  28th-order 


harmonics  are  expected  to  have  a much  smaller  effect.  Of  the  terms  with  y • 1 
those  with  q ■ 0 , 1 and  -1  are  likely  to  be  the  most  impqrtant,  since 
terms  with  q - ±2  have  an  extra  e factor,  and  e < 0.08  for  the  satellites 


analysed  here.  So,  for  the  present,  we  concentrate  on  the  terms  with 


^ c . 
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The  most  important  terms  in  (10)  are  likely  to  be  those  with  (y,q)  ■ (1,1)  and 
(1,  -1),  but  for  consistency  with  equation  (6)  we  also  give  explicitly  the  terms 
with  (y,q)  * (1,0)  . Equation  (10)  may  then  be  written  for  14th-order  resonance 
as 


£ ' ! (^'['(t)  sin  * * C ♦) 


14 
-1,0 
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»,  14,7^ 

14,14,6{514*  sin(*  + u)  " 514>2  cos  + w)} 
terms  in  jje  ^ *|q  - i (k  + q)e2  j gin  “ qw)]J  (II) 
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3 METHOD 


where  the  and  S(^  are  given  by  equations  (7)  to  (9). 


For  a number  of  years  we  have  been  receiving  US  Navy  orbital  data  on 
selected  low-drag  satellites  in  near-circular  orbits  close  to  14th-order  resonance. 
With  the  data  now  accumulated  on  two  of  these  satellites,  1965-16G  and  1971-120B, 
accurate  fitting  of  the  14th-order  theoretical  equations  is  possible.  Useful 
results  are  also  derived  from  the  US  Navy  data  on  1965-8IA,  a satellite  of  higher 
eccentricity,  which  passed  through  resonance  in  1970.  We  first  fit  the  inclina- 
tion and  eccentricity  data  separately,  using  the  computer  program  THROE and 

then  make  a combined  fitting  of  inclination  and  eccentricity,  using  the  S1MRES 
1 2 

program  . Our  fourth  satellite  is  Ariel  2 (1964-15A),  for  which  Gooding  has 

analysed  the  variations  in  inclination  and  eccentricity  after  14th-order 
12 

resonance 

In  determining  the  coefficients  of  14th  order  and  odd  degree  (degree  15, 

17,  19  ...),  we  have  utilized  data  from  six  other  satellites.  Of  these,  the 
most  important  are  three  low-drag  satellites  accurately  analysed  by  Wagner^: 
I963-26A;  196I-15G;  and  1971—1 20A  before  resonance.  The  other  three  are  high- 
drag  satellites,  two  analysed  by  Kloko2nik*^  and  one  by  Hiller*2,  for  which  the 
results  are  inevitably  of  poorer  accuracy,  but  still  useful  in  filling  gaps  in 
the  range  of  inclinations. 

In  determining  coefficients  of  14th  order  and  even  degree  (degree  14,  16, 

1 8 

18  ....),  the  four  primary  satellites  provide  the  main  data  base;  but  results 
from  the  high-drag  satellite  1971—1 06 A are  also  included,  and  some  subsidiary 
results  from  Wagner's  fittings  of  1963-26A  and  I961-15G  are  used. 
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THE  INDIVIDUAL  SATELLITES 


4. I The  four  primary  orbits 

4.1.1  Surcal  2,  1 965-1 6G 

Surcal  2 was  one  of  eight  small  satellites  launched  on  9 March  1965  into 
near-circular  orbits  having  inclinations  near  70°  and  heights  near  920  km.  Of 
the  eight  satellites,  1965-I6G  has  decayed  most  rapidly:  its  period  decreased 
from  103.5  minutes  initially  to  101.8  minutes  in  July  1977,  and  it  passed 
through  14th-order  resonance  on  13  November  1976. 

Values  of  inclination  from  US  Navy  orbits  from  January  1976  to  July  1977 

are  plotted  in  Fig  1,  after  removal  of  (1)  lunisolar  and  zonal  harmonic  perturba- 
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tions,  using  the  computer  program  PROD  with  1-day  integration  steps;  and  (2) 

the  effects  of  atmospheric  rotation  and  the  precession  of  the  pole,  using  THROE. 

The  unbroken  line  in  Fig  1 shows  the  theoretical  curve  fitted  by  THROE  with 

(y,q)  * (1,0),  while  the  broken  line  shows  the  fitting  of  inclination  and 

eccentricity  together  using  SIMRES,  with  (y,q)  * (1,0),  (1,1)  and  (1,-1).  The 

close  agreement  of  the  two  curves  is  extremely  satisfactory;  the  rms  scatter  of 

the  points  about  the  curve,  0.0012°,  is  slightly  smaller  than  the  expected 

accuracy  of  the  US  Navy  orbits.  The  values  of  the  lumped  coefficients  for 

(y,q)  * (1*0)  from  the  THROE  run  are  virtually  the  same  as  from  SIMRES.  Each 

2 

THROE  run  includes  a calculation  of  the  measure  of  fit  c (where  e is  the 
sum  of  squares  of  weighted  residuals,  divided  by  the  number  of  degrees  of  freedom). 
In  the  SIMRES  fitting,  i and  e were  weighted  on  the  basis  of  the  values  of 
c obtained  in  the  contributing  THROE  runs.  In  Fig  1 the  value  of  & increases 
from  -2.9  deg/day  at  the  start  to  +2.8  deg/day  at  the  end. 

The  values  of  eccentricity  from  US  Navy  orbits  are  plotted  in  Fig  2 after 
removal  of  zonal  harmonic  perturbations.  The  unbroken  line  shows  the  curve  fitted 
by  THROE  with  (y,q)  ■ (1,1)  and  (1,-1),  the  broken  line  shows  the  fitting  of  i 
and  e together  by  SIMRES,  with  (y,q)  » (1,0),  (1,1)  and  (1,-1),  and  the  dotted 
line  shows  the  THROE  solution  with  (y,q)  - (1,0)  (1,±1)  and  (2,±1).  The  agree- 
ment between  the  first  two  curves  is  good,  but  the  points  exhibit  a curious 
unfitted  oscillation  with  an  amplitude  of  order  0.0002  and  a period  of  about 
75  days.  This  cannot  be  due  to  lunisolar  perturbations,  which  have  been  calcu- 
lated and  amount  to  only  0.00001;  and  there  is  very  little  improvement  in  the 
fitting  on  taking  into  account  the  terms  (y.q)  - (1,±2)  or  (0,3).  But  the 
addition  of  the  ty,q)  ■ (2 , ± 1 ) terms  does  give  a better  fit:  the  measure  of  fit 
e decreases  from  1.10  to  0.94.  In  the  SIMRES  solution  the  rms  scatter  of  the 
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points  about  the  curve  is  0.00012,  which  is  rather  greater  than  the  expected 
accuracy  of  the  US  Navy  values. 

Because  of  these  deficiencies  in  fitting,  we  looked  with  some  suspicion  on 
the  values  of  the  lumped  coefficients  obtained  from  SIMRES  for  the  (y,q)  * (1,±1) 
terms.  This  suspicion  proved  well  founded:  in  the  solutions  for  individual 
coefficients  of  even  degree,  the  lumped  harmonics  did  not  conform  at  all  well 
with  those  obtained  from  the  other  satellites,  even  when  three  of  the  four  stan- 
dard deviations  were  increased  by  a factor  of  3.  Ue  therefore  looked  again  at 
the  THROE  runs  for  e and  found  that  the  run  with  (y,q)  * (1,0),  (l,±l)  and 
(2,±1)  gave  values  in  better  conformity  with  the  other  satellites.  The  values  of 
lumped  coefficients  for  q * ± 1 from  this  run  were  substituted  in  the  equations 

(without  any  modifications  in  sd)  and  gave  better  solutions  for  the  individual 

coefficients:  although  only  two  out  of  14  equations  were  altered,  the  standard 
deviations  of  the  solutions  were  reduced  by  I5Z  on  average.  Our  aim  is  to 
obtain  the  best  values  of  the  individual  coefficients,  so  this  ex  post  facto 
alteration  was  welcome.  This  fitting  is  presumably  better  because  the  q - ± 2 
(m  - 28)  terms  are  important  for  this  satellite,  and  should  not  be  allowed  to 
contaminate  the  q ■ ± 1 (m  « 14)  coefficients  which  we  are  trying  to  evaluate. 

Tables  I to  3 record  the  values  of  the  lumped  harmonics  as  used  in  the 
solutions  for  individual  coefficients,  for  I965-16G  and  all  the  other  satellites. 

The  values  for  1 965-1 6G  in  Table  1 are  from  the  SIMRES  fitting,  the  values  in 

Tables  2 and  3 from  the  THROE  fitting  of  e , as  already  mentioned.  These  Tables 

also  give  the  corresponding  values  of  the  coefficients  defined  in 

equation  (5),  and  similar  data  for  all  the  other  satellites  used. 

4.1.2  1971-120B  (Meteor  10  rocket?) 

The  Meteor  10  satellite  and  its  rocket  were  launched  on  29  December  1971 
into  a near-circular  orbit  at  a height  near  860  km,  at  an  inclination  of  81.25°. 

In  April  1976  three  fragments  appeared  in  orbit,  and  the  object  197I-120B,  which 
had  previously  been  approaching  14th-order  resonance,  reappeared  in  an  orbit  which 
had  just  passed  resonance.  Some  difference  of  opinion  exists  about  the  identity 
of  I97I-I20B,  but  this  is  irrelevant  in  our  analysis. 

Fig  3 shows  the  weekly  US  Navy  values  of  inclination  between  April  1976 
and  August  1977,  after  removal  of  the  usual  perturbations  (lunisolar,  odd  harmonic, 
etc).  The  unbroken  line  indicates  the  theoretical  curve  for  (y,q)  » (1,0) 
fitted  by  THROE.  The  rms  scatter  of  the  points  about  the  curve  is  0.0021.  The 
SIMRES  fitting  of  i and  e together,  with  equal  weighting  on  the  basis  of  the 


i 


Table 


Values  of  C.,  and  S..  from  the  (y*q)  = 0*0  terms  in  the  final  fittings,  as  used  in  the  solutions 


965-8 


12 


i 


THROE  runs,  is  shown  as  a broken  line  in  the  few  regions  where  it  differs  per- 
ceptibly from  the  unbroken  line.  The  value  of  ♦ increases  from  0.34  deg/day 
initially  to  0.58  deg/day  at  the  end.  The  fitting  of  inclination,  and  the 
agreement  between  the  two  curves,  is  excellent.  The  values  for  the  lumped 
coefficients  obtained  from  SIMRES  for  the  (y»q)  E 0,0)  terms  are 


9.0.' 

,0C14 


- 3.7  ± 1.4 


and 


9- 

10  S 


0,1 

14 


- 19.1  ± 1.2 


Wagner  obtained  values  of  these  lumped  harmonics  for  1971-120A  before  resonance. 
His  values  are  1.4  ± 1.2  and  -20.0  ± 1.2  . The  agreement  is  excellent,  con- 
sidering that  different  'wings'  of  the  resonance  are  being  analysed,  and  that 

neither  analysis  covers  the  exact  resonance.  (In  the  final  solution  the  sd  was 

-0,1 

doubled  on  both  values  of  C ^ , to  avoid  their  clashing  with  each  other.) 


The  values  of  eccentricity  for  1971—1 20B , after  removal  of  odd  zonal  harmonic 
perturbations,  are  plotted  in  Fig  4,  and  the  unbroken  line  shows  the  theoretical 
curve  for  (y,q)  = (1,1)  and  (1,-1)  as  fitted  by  THROE.  The  SIMRES  fitting  of  i 
and  e together  gives  an  almost  identical  result.  Fig  4 shows  that  the  fitting 
is  extraordinarily  good,  the  rms  scatter  of  the  points  about  the  curve  being 
0.000035,  which  is  smaller  than  the  likely  errors  in  the  values. 

4.1.3  OGO  2,  1965-81 A 

Launched  on  14  October  1965,  OGO  2 entered  an  eccentric  orbit  with  perigee 
height  420  km  and  apogee  1520  km,  at  an  inclination  of  87.4°.  Because  of  its 
low  perigee,  1 965-8 1A  suffered  much  greater  drag  than  the  two  satellites  previously 
discussed:  the  period,  initially  104.4  minutes,  steadily  decreased  under  the 
influence  of  air  drag,  and  the  orbit  passed  through  14th-order  resonance  on 
17  June  1970. 


In  equation  (6)  for  di/dt  , the  first  term  in  curly  brackets  (for  q ■ 0) 
is  usually  dominant.  But  1965-81A  happens  to  be  near  the  inclination  (86.2°) 
where  F , 7 is  zero;  also  the  eccentricity  is  quite  large,  0.064,  and  con- 
sequently  the  q = +1  terms  are  larger.  Numerically,  we  find 


cos  i)F 


15,14,7 


< 1 4)Fj  4, |4, 7^ 


0.2  . 


! 


I 
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-1,0  -1,0 

Also  it  happens  that  at  this  inclination  C,.  and  S,,  are  about  five  times 

-0,1  -0,1  ** 

larger  numerically  than  and  S14  . So  we  might  expect  that  the 

(Y,q)  - (1,0)  terms  in  (6)  would  be  negligible.  This  expectation  is  fulfilled. 
Fitting  i with  (y»q)  * (1,0)  alone  gives  poor  results;  a much  better  fitting 
is  obtained  with  (1,1)  and  (1,-1);  and  the  use  of  all  three  pairs  gives  no 
advantage. 

So  for  1 965-81 A,  both  the  inclination  and  eccentricity  are  best  fitted  by 
using  the  terms  (y,q)  * (1,1)  and  (1,-1)  only.  Both  orbital  parameters  therefore 
provide  values  of  the  same  lumped  coefficients,  and  this  satellite  provides  a 
severe  test  of  their  compatibility.  Table  4 gives  values  of  these  four  lumped 
coefficients  obtained  from  the  fitting  of  i alone  by  THROE,  e alone  by  THROE 
and  i and  e by  S IMRES.  The  agreement  is  as  good  as  can  be  expected  from  a 
fitting  of  four  coefficients  to  only  26  data  points. 

Figs  5 and  6 show  the  observational  values  and  the  fitted  curves.  The 
fittings  of  i in  Fig  5 are  quite  satisfactory,  but  Fig  6 gives  the  isq>ression 
that  the  fitting  of  e alone  is  rather  too  oscillatory,  and  that  the  SIMRES 
fitting  is  not  oscillatory  enough.  This  is  in  conformity  with  Table  4,  where  all 
the  values  from  the  SIMRES  solution  are  numerically  smaller  than  those  obtained 
from  fitting  e alone.  All-in-all,  the  fittings  and  the  lumped  coefficients  are 
good,  considering  the  small  number  of  observational  values.  In  Figs  5 and  6 the 
values  of  * run  from  -12.4  deg/day  at  the  beginning  to  +3.9  deg/day  at  the  end. 
The  rms  scatter  of  the  values  about  the  unbroken  curves  in  Figs  5 and  6 is  0.0012° 
for  i and  0.00007  for  e ; for  the  broken  curves  the  corresponding  values  are 
0.0013°  and  0.00012. 

For  1965-8IA  the  possibility  arises  that  the  terms  with  (y»q)  “ (I, ±2)  may 

be  important.  Their  inclusion  was  fruitless  with  i , but  did  improve  the  fitting 

for  e ; however,  too  many  constants  (nine  in  all)  were  being  determined  from 

too  few  points  (26),  so  it  was  no  surprise  to  find  that  nearly  all  the  lumped 

coefficients  were  undetermined  and  apparently  much  too  large.  To  assess  the 

- - 2 -1  - - -2  3 

effects  of  the  relevant  lumped  coefficients,  (C,S)j*  and  (C,S)j^*  , their 

values  were  computed  from  the  values  of  the  individual  coefficients  of  degree  15, 
17,  19  and  21  (Table  5),  and  these  values  were  used  with  THROE  to  calculate  the 
effect  of  the  (y,q)  ■ (1,±2)  terms  on  e . The  effect  was  found  to  be  very  small, 
the  maximum  change  in  e during  the  run  being  0.000012.  So  the  (y,q)  ■ 0,±2) 
terms  were  ignored,  and  we  adopted  the  SIMRES  solution  for  (y,q)  ■ (l,±l)  given 
in  Table  4.  In  the  final  solution  the  sd  of  S^’  was  doubled,  but  the  other 
015  values  fitted  well. 


1A 


Table  A 

Values  of  lumped  coefficients  obtained  from  I965-81A  by  fitting  i alone, 

e alone,  or  i and  e together 


From  i alone 
THROE 

From  e alone 
THROE 

From  e and  i 
SIMRES 

9J’° 

,0C1A 

-55  ± 7 

-152  ± 18 

-60  ± 7 

.1,0 

S1A 

-32  ± 9 

-92  t 26 

-35  ± 11 

--I.2 

C1A 

-A3  ± 7 

-63  ± 18 

-52  ± 8 

-->.2 
S,  A 

9 ± 17 

29A  ± 12 

12  ± 17 

A. I .A  Ariel  2,  196A-15A 

Ariel  2,  the  second  Anglo-US  satellite,  was  launched  on  27  March  196A  into 

an  orbit  with  inclination  51.6A°  with  period  101.3  minutes.  Initially,  the 

heights  of  perigee  and  apogee  were  290  and  1360  km  respectively,  and  the  eccen- 

20 

tricity  was  0.075.  Using  Minitrack  observations,  Gooding  computed  the  orbit 

12 

at  210  epochs  during  the  first  year  after  launch,  and  has  now  determined  values 
of  lumped  lAth-order  harmonics  by  analysing  the  variations  in  inclination  and 
eccentricity;  the  orbit  was  already  post-resonant  at  launch.  Because  of  the 
relatively  high  eccentricity  and  the  multitude  of  values,  the  best  fitting  is 
obtained  with  eight  pairs  of  values  of  (y,q)  - the  basic  trio  and  five  other 
paira.  The  values  of  inclination  and  eccentricity,  and  the  curves  fitted  by 
SIMRES,  are  shown  in  Figs  7 and  8.  Although  the  fittings  do  not  look  so  convin- 
cing as  those  in  Figa  1 to  6,  there  are  far  more  values  available  for  Ariel  2,  and 

the  variations  in  the  orbital  parameters  extend  much  further  away  from  resonance, 

2 1 

because  of  the  high  eccentricity,  as  also  happens  for  example  with  Vanguard  3. 

For  Ariel  2,  the  Q coefficients  do  not  fall  off  so  rapidly  as  for  the 

three  previous  satellites  (see  Tables  I to  3),  and  the  neglect  of  harmoqics  of 

degree  greater  than  23  may  be  a source  of  error.  The  effect  of  these  higher 

harmonics  was  calculated  using  the  values  in  Tables  5 and  7,  and  the  only  sd 

-0,1 

needing  to  be  increased  for  this  reason  was  that  of  C,,  , which  was  increased 
-9 

from  23  to  30  x 10  . The  numerical  values  of  the  lumped  harmonics  as  used  in 

the  solutions  are  given  in  Tables  I to  3. 


015 
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4.2  The  orbits  analysed  by  Wagner 

In  a thorough  assessment  of  the  accuracy  of  the  Goddard  Earth  Models, 

Wagner"*  has  analysed  the  l4th-order  resonant  variations  in  inclination  (but  not 
eccentricity)  of  three  satellites  relevant  to  our  work.  The  first  of  these  three 
orbits  is  that  of  I963-26A,  which  passed  through  l4th-order  resonance  in  March 
1968,  with  inclination  49.7°  and  eccentricity  0.057.  The  second  orbit  is  that 
of  1961-15G,  a fragment  from  the  exploded  rocket  of  Transit  4A;  this  fragment 
passed  through  l4th-order  resonance  in  September  1971,  with  inclination  66.6°  and 
eccentricity  0.017.  The  third  satellite  is  1971-120A,  Meteor  10,  which  Wagner 
• analysed  for  the  years  1972-1975,  some  years  before  it  was  due  to  reach  resonance. 

As  already  mentioned,  Wagner's  values  are  close  to  those  we  have  obtained  from 
1971-I20B  after  resonance,  and  the  two  satellites,  which  are  in  nearly  identical 
orbits,  should  together  define  a reliable  mean  value. 

Wagner's  values  for  the  q*  0 coefficients  from  these  three  satellites 

-0,1 

are  given  in  Table  1,  the  standard  deviation  for  for  1971-120A  being 

doubled,  as  mentioned  in  section  4.1.2. 

A fourth  orbit  analysed  by  Wagner,  or  rather  a set  of  orbits  of  several 
fragments  from  the  1965-82  launch  at  inclination  32°,  has  been  used  to  strengthen 
our  8-coefficient  solution  for  odd-degree  coefficients.  For  1965-82  the  values 
of  the  are  largest  for  l ■ 23  to  29,  and  the  geopotential  coefficients  of 

these  degrees  will  make  a dominant  contribution  to  the  lumped  coefficient.  So 
1965-82  is  unsuitable  unless  at  least  eight  coefficients  are  being  evaluated. 
Numerical  data  appear  at  the  end  of  Table  1. 

For  I963-26A  and  1961-15G,  Wagner  also  obtains  values  for  the  q - ±1 
coefficients.  These  are  given  in  Tables  2 and  3,  but  the  standard  deviations 
obtained  by  Wagner^  have  been  doubled,  because  his  values  are  derived  from  the 
, variation  in  inclination  only,  and  are  not  combined  inclination-eccentricity 

fittings,  as  our  q ■ i 1 values  are. 

4.3  The  four  high-drag  orbits 

KlokoJnik16  has  analysed  the  14th-order  resonant  variations  in  inclination 
(but  not  eccentricity)  for  two  satellites  with  perigee  heights  near  200  km, 
namely  Intercosmos  9,  1973-22A,  inclination  48.4°,  and  Intercosmos  10,  1973-82A, 
inclination  74°.  Because  of  the  high  drag,  the  values  of  the  lumped  harmonics 
are  rather  inaccurate,  and  in  the  solution  it  was  necessary  to  increase  the  sd 
: 0 1 5 of  1973-22A  by  a factor  of  2 and  the  sd  of  I973-82A  by  a factor  of  5. 


A 
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Hiller  has  recently  analysed  the  14th-order  resonant  variation  of  inclina- 
tion and  eccentricity  for  China  2 rocket,  1971-18B,  at  inclination  69.9°,  but 

again  the  orbit  did  not  remain  near  resonance  for  long  enough  to  allow  an  accurate 

— 0 , 1 

analysis,  and  it  was  necessary  to  increase  the  sd  of  by  a factor  of  2. 

Despite  their  limitations,  the  results  from  these  three  orbits  helped  to 
fill  gaps  in  the  coverage  of  inclination,  in  the  odd-degree  solutions.  The 
values  of  the  lumped  coefficients  from  the  three  orbits  are  given  in  Table  1. 

The  fourth  high-drag  orbit  used  is  that  of  Cosmos  462,  1971-106A,  with 

1 8 

perigee  height  240  km,  which  has  been  analysed  by  Walker  . This  orbit  has  some 
similarities  with  1965-81 A.  The  q ■ 0 terms  are  found  to  have  little  effect,  and 
the  analysis  yields  values  of  the  q - +1  coefficients,  which  are  given  in  Tables  2 
and  3.  In  the  analysis  of  1971-106A,  the  perturbations  in  i and  e due  to 
q ■ ±2  terms  were  calculated,  using  lumped  coefficients  evaluated  from  the  solu- 
tions of  Table  5,  and  the  values  of  i and  e were  then  cleared  of  these  per- 
turbations. This  is  a refinement  which  should  always  be  used  in  resonance  analyses 
if  reliable  lumped  harmonics  can  be  pre-calculated. 

4.4  Data  not  used 

Many  satellites  move  in  orbits  which  are  close  to  14th-order  resonance  but 

do  not  experience  exact  resonance.  These  orbits  in  'shallow  resonance'  (as  it 

is  often  called)  have  been  extensively  studied  over  the  years  and  utilized  in 

22 

the  comprehensive  geoid  models.  Also  Reigber  and  Balmino  have  obtained  solu- 
tions for  the  individual  coefficients  of  14th  order  up  to  degree  30,  from  analysis 
of  a number  of  shallow  resonant  orbits.  Their  results  have  been  incorporated  in 
the  GRIM  2 model;  so  we  thought  it  best  not  to  mix  their  values  with  ours,  but  to 
obtain  an  independent  solution  for  comparison. 

5 THE  SOLUTIONS 

5.1  Solutions  for  odd-degree  coefficients  (*.  » 15,  17.  19  ...) 


The  fittings  of  the  (y,q)  ■ (1,0)  terms  give  us  nine  values  of  the 
lumped  coefficients  with  [q,k]  - [0,1],  listed  in  Table  I.  Inserting  these 
values  in  equations  (5),  we  obtain  nine  equations  for  14th-order  C-coeff icients 
of  odd  degree, 

- 0.1-  0 1-  .0,1 

u * *[|7  u + 9ia  u + ••••  “ (15 


''15,14  T XI7  ''l?, 14  T yl9  ^19, 14  T **14  ' 

-0,1 

with  nine  similar  equations  for  the  coefficients  S ^ . The  numerical  values 

of  the  Qg  terms  and  lumped  coefficients  for  each  satellite  are  given  in  the 
first  nine  rows  of  Table  I. 
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We  follow  the  method  which  proved  successful  in  solving  for  l5th-order 
coefficients.  To  the  nine  equations  of  the  form  (12)  we  add  constraint  equations 


of  the  form 


0 ± I0-5/*2 


and  then  solve  by  least  squares  for  3,  4,  5 . ...r  coefficients  from  12,  13,  14  .... 
14  ....(9  + r)  equations.  The  value  of  r is  chosen  empirically  after  examina- 
tion of  the  solutions,  to  give  what  is  believed  to  be  the  best  set  of  values  for 
the  individual  coefficients.  The  quantity  10  /*.  derives  from  "Raula's  rule 

of  thumb"  for  the  order  of  magnitude  of  a coefficient  of  degree  l . Now  that 
individual  coefficients  are  being  evaluated  with  better  accuracy,  it  is  becoming 

apparent  that,  at  least  for  12  < £.  < 30  , a better  average  value  for  the 

—5  2 —5  2 

coefficients  would  be  about  0.7  * 10  /£  , with  10  /l  itself  regarded  rather 

as  an  upper  limit.  Ue  found  that  the  standard  deviations  of  our  solutions  were 
slightly  improved  by  making  this  30Z  reduction  in  the  right-hand  side  of 
equation  (13);  but  the  values  of  the  coefficients  were  little  altered,  so  it 
seemed  better  to  allow  the  more  relaxed  constraint,  as  given  in  (13). 

In  the  course  of  the  solutions  of  equations  (12)  and  (13),  some  values  of 

lumped  harmonics  were  found  to  be  persistently  ill-fitting.  We  tried  doubling 
their  assumed  standard  deviations,  and  if  the  solution  was  significantly  improved, 
the  doubled  sd  was  used*.  The  changes  made  have  already  been  mentioned;  the 
values  given  in  Table  I incorporate  the  increased  sd  and  are  those  used  in  the 
solutions. 

When  the  nine  equations  itemized  in  Table  I,  together  with  the  constraint 
equations  (13),  were  solved  by  least  squares  for  3,  4,  5 ....9  coefficients,  the 
values  of  e , the  measure  of  fit,  were  2.13,  1.08,  1.07,  1.07,  1.06,  1.06,  1.06 

for  the  C-coefficients,  and  1.23,  0.59,  0.59,  0.58,  0.58,  0.58,  0.58,  for  the 

. . 2 

S-coefficients.  (As  before,  e is  defined  as  the  sum  of  the  squares  of  the 
weighted  residuals  divided  by  the  number  of  degrees  of  freedom,  which  here  is 
always  9.)  Since  e fails  to  decrease  when  the  number  of  coefficients  evaluated 
goes  beyond  4,  the  4-coefficient  solutions,  for  «.  - 15,  17,  19  and  21,  are  recom- 
mended as  best,  and  are  underlined  in  Table  5. 

* Purists  sometimes  query  the  legitimacy  of  this  procedure.  Our  aim  is  to  derive 
the  best  possible  values  of  the  geopotential  coefficients:  this  aim  is  not 
helped  by  giving  too  high  a weight  to  inaccurate  data.  The  inaccurate  values 
could  be  omitted  altogether;  but  this  would  imply  a judgment  that  zero  weighting 
is  optimal,  which  is  unlikely  to  be  true  unless  the  data  are  entirely  spurious. 
Our  procedure  is  an  attempt  to  optimize  the  weighting  empirically. 
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The  8-coefficient  solutions  are  also  given  in  Table  5,  however,  for  two 
reasons:  (a)  to  show  how  the  values  of  the  first  four  coefficients  are  affected 
by  inclusion  of  higher  harmonics;  and  (b)  to  provide  a comparison  with  GEM  10  and 
GRIM  2,  which  go  to  degree  29.  When  solving  for  eight  coefficients,  we  include 
a further  equation,  that  from  the  1965-82  satellites  at  32°  inclination,  given 
in  the  last  row  of  Table  1 . 

Table  5 

Solutions  for  geopotential  harmonics  of  order  14  and  odd  degree 


OIS 
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The  weighted  residual,  that  is 

(observational  value  minus  computed  value)  * (standard  deviation), 
for  each  satellite  in  the  4-coeff icient  solution  is  given  in  Table  6.  The  resi- 
duals for  the  8-coefficient  solution  are  very  similar.  It  is  obvious  from 
Table  6 that  the  formal  standard  deviations  of  the  solutions  for  the  C-coefficients 
could  be  reduced  by  omitting  the  three  high-drag  satellites.  When  the  equations 
are  solved  with  these  three  satellites  omitted,  it  is  found  that  the  values  of  the 
C-coefficients  remain  within  { sd  of  those  in  Table  5,  but  the  standard  deviations 
are  reduced  by  20Z  on  average.  Although  this  solution  is  formally  more  accurate, 
it  is  probably  less  reliable  because  of  the  wider  gaps  in  the  coverage  of  inclina- 
tion, so  we  prefer  the  solution  given  in  Table  5. 


Equation  (6)  shows  that  a truer  measure  of  the  strength  of  the  resonance 
effects  is  given  not  by  the  lumped  coefficients  alone  but  by  these  coefficients 
multiplied  by  F , 7 , and  values  of  the  lumped  coefficients  multiplied  in 
this  way  are  given  in  the  last  two  columns  of  Table  1 (which  also  shows  that  the 
most  accurate  results  are  those  from  1961-1SG,  1965-16G  and  1971-120). 

The  values  of  the  lumped  coefficients  multiplied  by  F^  ^ ^ are  also 

plotted  in  Fig  9,  with  curves  showing  the  values  given  by  the  4-coefficient 

solutions.  It  is  seen  that  the  solutions  are  entirely  satisfactory,  in  the 

sense  that  there  are  no  strong  oscillations  in  the  gaps  between  the  data  points. 

The  points  indicated  by  the  triangles  in  Fig  9 are  at  an  inclination  of  86.2°, 

where  F,r  ,,  , * 0 : a satellite  in  a circular  orbit  at  this  inclination  would 
15,14,7 

suffer  no  perturbation  in  inclination  as  it  passes  14th-order  resonance  (except 
possibly  from  28th-order  coefficients). 

5.2  Solutions  for  even-degree  coefficients  (£  = 14,  16,  18  ...) 

For  the  orbits  of  the  four  basic  satellites  and  1971-106A,  the  (y,q)  * ( 1 ,± 1 ) 
terms  in  the  fittings  of  the  variations  in  eccentricity  and  inclination  have  given 
two  equations  each  for  the  !4th-order  coefficients  of  even  degree,  of  the  form 


^14,14 

^14,14 


. nl ,0-  _1 ,0- 

q!6  C16, 14  Q18  C 1 8 , 14 


+ 0- ^ * 2c  ♦ o' ^ 1 2c  + 
qi6  C 16, 14  qi 8 C18, 14  •••• 


C 


1,0 

14 


04) 

(15) 
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and  similarly  for  S , on  replacing  C by  S throughout.  The  numerical  values 
of  the  Qg  and  the  lumped  harmonics  on  the  right-hand  side,  for  each  of  these 
satellites  and  for  the  two  satellites  of  Wagner,  are  given  in  Table  2,  for 
equation  (14),  and  Table  3,  for  equation  (15).  To  these  14  equations  for 

^ (and  14  for  ^)  we  add  the  constraint  equations  (13),  and  solve  for 
r coefficients  from  (14  + r)  equations,  with  r - 3,  4 9 . 


With  14  equations  to  be  fitted  simultaneously,  there  is  little  scope  for 

the  values  of  the  coefficients  to  adjust  themselves  to  fit  non-conforming  lumped 

values.  So  it  was  no  surprise  to  find  a number  of  obstinately  ill-fitting 

lumped  coefficients,  for  which  it  was  advantageous  to  increase  the  sd  . They 

_~1 , 2 


were  as  follows.  The  sd  of  C..  was  doubled  for  196I-I5G  and  1 971-1 20B, 

14  _1,0 

and  multiplied  by  10  for  1971-106A.  The  sd  of  S..  was  multiplied  by  5 for 

-12 

1971-106A.  The  sd  of  3,  ’ was  doubled  for  1963-26A  and  1965-81 A.  The 

14 

values  given  in  Tables  2 and  3 are  those  used  in  the  solutions,  after  these 
increases  in  sd  have  been  made. 


The  14  equations  of  type  (14)  or  (15),  itemized  in  Tables  2 and  3,  plus 
the  constraint  equations  (13),  were  solved  by  least  squares  for  3,  4 ....9 
coefficients.  The  values  of  e , the  measure  of  fit,  were  1.03,  1.01,  0.91,  0.89, 
0.89,  0.89,  0.89  for  the  C-coef ficients  and  0.92,  0.74,  0.69,  0.68,  0.68,  0.68, 
0.68  for  the  S-coeff icients.  For  both  C and  S , therefore,  e decreases 
significantly  - by  102  and  11  respectively  - on  going  from  the  4-  to  the 
5-coefficient  solution,  but  does  not  decrease  significantly  thereafter.  So  the 
5-coefficient  solutions,  for  £ = 14,  16,  18,  20  and  22,  are  recommended  as  the 
best,  and  are  given  underlined  in  Table  7.  The  9-coefficient  solution  is  also 
given  in  Table  7,  however,  to  show  how  little  the  solution  is  affected  by 
including  higher-degree  coefficients,  and  to  provide  a comparison  with  GRIM  2. 

Equations  (6)  and  (31)  show  that  a bejtt^r  measure  of  the^e^fect  of  the 
lumped  coefficients  is  given  by  ^ and  F ^ ^ than  by  the 

lumped  coefficients  alone.  The  values  of  these  quantities  (S  as  well  as  C)  are 
listed  in  the  last  two  columns  of  Tables  2 and  3,  and  the  values  are  plotted  in 
Figs  10  and  11,  with  the  curves  given  by  the  5-coefficient  solution.  From  Fig  11 
alone  it  may  seem  at  first  sight  that  the  curves  are  making  rather  a poor  job  of 
fitting  the  points;  but  the  fit  is  of  course  simultaneous  for  the  C (or  S) 
values  in  both  Fig  10  and  Fig  11,  so  that  the  form  of  the  curve  is  a compromise 
between  their  conflicting  claims. 


Table  7 


Solutions  for  geopotential  harmonics  of  order  14  and  even  degree 


Table  8 

Weighted  residuals  for  each  satellite  in  the 
5-coefficient  solution  of  Table  7 

Residuals  for 


Satellite 

_1,0 

CI4 

-1,0 

S.4 

--U2 

CI4 

c"1’2 

S|4 

1 963-26A 

0.49 

-0.29 

0.64 

-0.95 

1 964-1 5A 

-0.72 

0.54 

-0.26 

-0.15 

1971-106A 

0.04 

-1.17 

-1.14 

0.71 

1 96 1—1 5G 

-0.07 

0.18 

-1.43 

-0.07 

1965-I6G 

0.68 

0.30 

-0.76 

-0.51 

I971-120B 

0.45 

0.01 

-1.75 

-0.08 

1965-81 A 

0.13 

-0.31 

0.77 

1.20 

6 COMPARISONS 

6.1  Odd-degree  coefficients 

Comparison  between  the  4-coefficient  and  the  8-coefficient  solutions  in 
Table  5 shows  that  the  solutions  are  extremely  stable.  The  values  in  the 
8-coefficient  solution  differ  from  the  corresponding  values  in  the  4-coefficient 
solution  by  less  than  the  standard  deviation  of  the  latter  (except  for  S|5  ^ , 
where  the  sd  is  unusually  small).  This  stability  is  a very  welcome  feature  of 
the  solutions. 
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Comparison  of  the  4-  and  8-coefficient  solutions  with  GEM  10  shows  good 
agreement  among  the  C coefficients  (within  1.5  sd);  but  our  S coefficients, 
which  have  very  low  standard  deviations,  do  not  agree  so  well  with  GEM  10  (except 
for  l m 15,  where  the  difference  is  only  3Z) . 

Comparison  with  GRIM  2 shows  poor  agreement  with  C coefficients  and  better 
agreement  than  GEM  10  for  the  S coefficients. 

Comparison  with  SSE  IV. 3 shows  quite  good  agreement  with  our  solutions  for 
the  C coefficients,  but  considerable  divergences  for  the  S coefficients. 

Our  value  and  those  from  the  three  comprehensive  geoid  solutions  are 
largely  independent,  but  not  entirely  so,  since  Wagner's  results  have  some 
influence  on  GEM  10,  as  well  as  on  our  values. 

The  GEM  8 model  agrees  with  our  values  just  about  as  well  as  GEM  10.  How- 
ever, no  detailed  comparisons  are  made,  because  we  feel  that  the  strongly 

oscillatory  high-degree  terms  in  GEM  8 for  the  C coefficients  (the  values  are 

-9 

-27,  +26  and  -41  * 10  for  i **  25,  27  and  29)  are  more  likely  to  be  an  artefact 
of  the  solution  than  realistic  values.  This  oscillatory  tendency  may  also  affect 
GEM  10  and  SSE  IV. 3,  though  to  a lesser  extent. 

Although  the  three  comprehensive  models  do  not  always  agree  either  with 
each  other  or  with  our  values,  they  agree  much  better  with  our  values  than 
happened  with  the  15th-order  coefficients,  no  doubt  because  several  shallow 
resonant  14th-order  orbits  are  included  in  the  comprehensive  solutions.  The 
agreement  is  best  for  t ■ 15,  17  and  19,  with  all  four  solutions  indicating: 

a low  positive  value  (around  5*)  for  ^ and  a value  close  to  -24  for 

c . 

15,14  ’ 

a value  near  -15  for  C,_  ..  , and  near  +15  for  S,_  ..  ; 

17,14  17,14 

for  C)9  ^ , the  verdict  is  'small  and  probably  negative',  and  for 

S,9  ^ 'negative  and  larger'  (near  -10). 

Ft  higher  values  of  l the  agreement  is  not  so  good,  and  the  near-unanimity 
over  and  ^ may  be  a chance  effect. 

The  comparison  may  perhaps  fairly  be  summarized  by  saying  that  our  solutions 
and  those  from  the  comprehensive  models  are  recognizably  similar,  but  further  work 
on  both  methods  is  needed  to  achieve  better  agreement. 
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6.2  Even-degree  coefficients 

Comparison  of  our  5-coefficient  and  9-coefficient  solutions  in  Table  7 
shows  that,  as  for  the  odd-degree  solutions,  the  values  are  extremely  stable. 

The  values  in  the  9-coefficient  solution  differ  from  the  corresponding  values 
in  the  5-coefficient  solution  by  less  than  half  the  standard  deviation  of  the 
latter.  The  chief  imperfection  is  the  rather  large  standard  deviations,  which 
are  two  or  three  times  greater  than  for  the  odd-degree  coefficients.  This 
reflects  the  fact  that  the  most  easily  measured  orbital  change  at  resonance,  that 
in  the  inclination,  usually  depends  primarily  on  the  odd-degree  coefficients. 

Comparison  of  our  9-coefficient  solutions  with  the  comprehensive  models 
shows  considerable  differences  in  the  C terms  for  l * 14  and  18.  For  l * 14 
the  values*  from  GEM  10,  GRIM  2 and  SSE  IV. 3 are  between  -50  and  -70,  whereas 
our  solution  gives  -39  ± 3.  For  1=  16,  the  three  comprehensive  models  differ, 
but  their  average  is  -23,  in  agreement  with  our  values.  For  l = 18  our  value 
is  -22 ±6,  while  the  three  comprehensive  models  give  values  between  -10  and  +1. 

For  the  S coefficients  the  agreement  is  much  better,  and  GEM  10  agrees  with  our 
values  to  within  1.4  sd  for  all  values  of  i (up  to  28).  For  l = 14  our 
5-coefficient  S solutions  give  -8 ±2,  while  the  other  three  range  between  -5 
and  -1.  For  l m \b  we  have  -36 ± 4,  while  the  other  three  range  between  -41  and 
-33.  For  18  the  agreement  is  not  so  good.  For  l ■ 20  the  GRIM  2 value  of 

-44.5  seems  unlikely,  but  the  other  two  are  fairly  close  to  ours.  For  X,  = 22  all 
three  are  within  2 sd  of  ours. 

For  even-degree  coefficients,  our  values  and  those  from  the  comprehensive 
models  are  almost  completely  independent,  since  Wagner's  limited  results  for 
(y,q)  “ (1**1)  have  very  little  influence  on  our  solutions. 

The  conclusions  from  Table  7 are  (a)  that  the  agreement  is  already  good 
for  , but  there  are  discrepancies  for  C^  ^ , and  (b)  that  the  standard 

deviations  of  our  solutions  need  to  be  improved  by  analysis  of  further  low-drag 
orbits. 

7 CONCLUSIONS 

We  have  for  the  first  time  determined  the  values  of  the  coefficients  of 
individual  l4th-order  harmonics  in  the  geopotential  from  analysis  of  changes  in 
the  inclination  and  eccentricity  of  satellites  with  orbits  which  passed  through 
l4th-order  resonance.  Results  from  the  eleven  orbits  are  used.  Although  values 
up  to  degree  30  are  obtained,  the  recommended  set  of  values  runs  from  degree  14 
to  22  only,  and  is  as  follows: 
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* In  this  paragraph  we  drop  the  factor  10~^ 
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-38.5  ± 2.9 
4.5  ± 1.1 
-22.3  ± 3.6 
-15.0  ± 2.6 
-24.0  ± 4.9 
-1.6  ± 2.8 
8.8  ± 5.8 
18.2  ± 3.6 
-14.5  ± 8.1 


-7.8  ± 2.2 
-23.8  ± 0.3 
-36.0  ± 3.8 
16.8  ± 1.2 
-3.2  ± 3.7 
-7.6  ± 1.0 
-15.4  ± 4.6 
-10.6  ± 1.9 
9.9  ± 6.4 


Tables  5 and  7 show  that  these  values  are  not  appreciably  disturbed  when  higher- 
degree  coefficients  are  included  in  the  evaluations.  The  odd-degree  coefficients, 
derived  primarily  from  changes  in  orbital  inclination,  are  more  accurate  than  the 
even-degree  coefficients,  which  are  determined  primarily  from  changes  in 
eccentricity. 

Values  obtained  in  this  way  from  resonance  effects  provide  an  independent 
test  of  the  accuracy  of  the  comprehensive  geoid  models  derived  in  recent  years. 
Detailed  comparisons  with  three  of  these  models,  GEM  10,  GRIM  2 and  SSE  IV. 3,  show 
that  their  values  of  l4th-order  harmonic  coefficients  are  generally  similar  to 
ours:  see  Tables  5 and  7.  The  best  agreement  - within  1$  sd  for  the  even-degree 
S coefficients  and  the  odd-degree  C coefficients  - is  with  GEM  10.  However, 

there  are  also  significant  discrepancies,  and  further  work  on  both  methods  is 
needed  to  achieve  better  agreement. 
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Fig  1 1966-16G:  variation  of  inclination  at  14th-order 


rom  US  Navy  orbits 
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Fig  2 1965-16G:  variation  of  a at  14th-order  resonance 


TR7M1B 


Fig  4 1971-120B:  variation  of  a near  14th -order 


from  US  Navy  orbits 
THROE  fitting  of  « alone  (1,1)  (1,-1 ) 


Ths  Barth's  gravitational  potential  ia  now  usually  expressed  in  tana  of  a 
double  series  of  tessera!  harmonic*  with  several  hundred  terns,  up  to  order  and 
degree  at  least  20.  the  harmonics  of  order  14  can  be  evaluated  by  analysing  changes 
ia  satellite  orbits  which  experience  1 4th-ordar  resonance,  when  the  track  over  the 
Berth  repeats  after  14  revolutions. 

H*  ■ A: 

■*'“*»»  tftis  Report  a)  describefew^first  evaluation  of  individual  14th-order 
coefficients  in  the  geopotential  frea  analysis  of  the  variations  in  inclination  end 
eccentricity  of  satellite  orbits  passing  through  !4th~order  resonance  under  the 
action  of  air  drag.^dteing  results  from  eleven  satellites  r we  iing£veluc8  for 
eighteen  normalised  coefficients  of  harmonics  of  order  14  and  degree  14,  15  ....22. 
These  values  provide  a test  of  the  accuracy  of  the  Uth-order  coefficients  in 
comprehensive  gaoid  models.  Detailed  comparisons  with  three  recent  modele  ere  made, 

coefficients  end  discrepancies  on  others*.. 


As  far  as  poasfele  this  page  dtotrid  contain  oaiy 
above  must  be  muted  to  indicate  the  daaatfteai 


1.  DR1C  Reference 
(to  bt  added  by  DRIC) 


\wvr~rrr-~m 

RAE  TR  78015 

M/A 

4.  Report  Security  Clamification/Maiking 
UNCLASSIFIED 


5.  DRIC  Code  for  Originator  4.  Oignstor  (Cotporaaa  Author)  Name  and  Location 

850100  Royal  Aircraft  Establishment,  Farnborough,  Hants*  UR 


5 a.  Sponsoring  Agency’s  Code 
N/A 


7.  Ihle  ||^ 

Evaluation  of  14th-order  harmonics  in  the  geopotential 


7a.  (For  Translations)  Title  in  Foreign  Language 


7b.  (For  Conference  Papers)  Title,  Place  and  Date  of  Conference 


8.  Author  1. 


i,  Initials  1 9a.  Author  2 


9b.  Authors  3, 4. 


King-Hsla*  D.G.  ( Walker,  Doreen  M.C.  | Gooding,  R.H. 


13.  Project 


11.  Contract  Number  12.  Period 

H/A  M/A 


15.  DMributioa  statement 

(a)  Controlled  by-  MOD(PE)  D.  Space 

(b)  Special  lioftattoaa  (if  any)  - 


14.  Descriptors  (Keywords)  (Descriptors  narked  * are  (elected  front  TEST) 

Gravitational  field.  Geopotential.  Satellite  orbits.  Resonance. 
Orbital  analysis.  . .1  , 5 i V 


10.  Date  Reft. 

February  33  „ 

1978  I 1 a 


14.  Other  Reference  Not. 
Spsce  545 


